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ON PANEL FLUTTER AND DIVERGENCE OF INFINITELY LONG UNSTIFFENED AND 
RING-STIFFENED THIN-WALLED CIRCULAR CYLINDERS 1 


By Robert W. Leonard and John M. Hedgepeth 


SUMMARY 


A ’preliminary theoretical investigation of the panel flutter and 
divergence of infinitely long, unstiffened and ring-stiffened 
thin-walled circular cylinders is described. Linearized un- 
steady potential-flow theory is utilized in conjunction with 
Donnell’s cylinder theory to obtain equilibrium equations for 
panel flutter. Where necessary, a simplified version of 
Flugge’s cylinder theory is used to obtain greater accuracy. 
By applying Nyquist diagram techniques, analytical criteria 
for the location of stability boundaries are derived. A limited 
number of computed results are presented. 

INTRODUCTION 


Although considerable effort has been expended in study- 
ing the flutter and divergence of thin, flat panels exposed to 
an airstream (see, for example, refs. l to 9), little is known 
of the importance of similar aeroelastic phenomena in the 
design of thin-walled cylindrical missile bodies or of other 
aircraft components where thin, curved panels are used. 
The purpose of this report is to describe a preliminary 
theoretical investigation of the aeroelastic stability of such 
configurations. Analytical criteria for the determination of 
panel-flutter and panel-divergence boundaries for infinitely 
long, unstiffened and ring-stiffened thin-walled circular 
cylinders are presented along with a limited number of 
computed results. 

SYMBOLS 
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Ci 
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e 

F,h\ 


distance botween ring stiffeners 
amplitude of mth term in expansion for lateral 
motion of ring-stiffened cylinder 
dummy variable 
speed of sound in air 
speed of sound in fluid inside cylinder 

speed of sound in cylinder material, -* — 

V Pi 


plate flexural stiffness per unit length, ■ 

Li n ) 

Young’s modulus 

base of natural system of logarithms 
outside and inside air-force functions defined by 
equations (5) and (6), respectively 
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Po,Pi 

Po,Pi 


R 

Ke 
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outside and inside air-force functions for vi- 
brating ring (see eqs. (34) and (35), respec- 
tively) 

Hankel functions of first and second kind, 
respectively, of order n 
thickness of cylinder wall 
imaginary part 

modified Bessel function of first kind of order n 

Bessel function of first kind of order n 
integer 

modified Bessel function of second kind of order n 
dimensionless frequency of harmonic vibration, 

— for unstiffened cylinder and — for ring- 

stiffened cylinder 

resonant frequencies of internal air-force func- 
tion for unstiffened cylinder (see eq. (22)) 
function defined by equation (8) 
function defined by equation (9) 
outward lift force on cylinder wall 
amplitude of harmonically varying outward 
lift force for unstiffened cylinder 
amplitude of mth term in expansion for outward 
lift force for stiffened cylinder 
Mach number of flow along vibrating cylinder 
Mach number of external or internal flow along 
stationary cylinder 

integer, number of longitudinal half-waves in 
each bay of ring-stiffened cylinder 
function defined by equation (44) 
integer, number of full waves around circum- 
ference of cylinder 

forces exerted on cylinder wall by ring stiffeners 
amplitudes of harmonically varying reaction 
forces exerted by ring stiffeners 
radius of cylinder 
real part 
radial coordinate 
time 

lateral deflection of cylinder wall, positive 
outward 


• Supersedes N AOA. Technical Note 363S by Robert W. Leonard and John M. Hedgepeth, 1056. 
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w amplitude of lateral motion of cylinder wall for 

unstiffened cylinder 

x longitudinal coordinate for vibrating cylinder 

Xi longitudinal coordinate for stationary cylinder 

z argument of Bessel functions for supersonic 

relative flow, ( r|V£ 2 — 1 for unstiffened cyl- 
inder and \m\a^i?—\ for stiffened cylinder 
a aspect-ratio parameter, Rvja 

function defined after equation (39) 

8(x) Dirac delta function, 5(z)=0 for 

f 5(z)dz=l 

c/ — co 

e damping coefficient 

£ argument of Bess el fimctions for subsonic relative 

flow, |7|Vl — £ 2 for unstiffened cylinder and 
|m]aVl — for stiffened cylinder 
0 angular circumferential coordinate 

n Poisson’s ratio 

v wave number of longitudinal waves in unstiffened 

cylinder 

v dimensionless wave number, Rv 

£ dummy variable 

p mass density of air 

Pi mass density of fluid inside cylinder 

p, mass density of cylinder material 

c rfl,<r x midplane stresses in circumferential and longi- 

__ tudinal directions, positive in tension 

<re,<r x dimensionless midplane stresses, c rg/E and <r s /E, 

respectively 

0 parameter defined after equation (8) 

Q m parameter defined after equation (39) 

« frequency of harmonic vibration 

V* operator, 

V -4 inverse of operator v 4 

Subscripts: 

cr critical value 

min minimum value 

max maximum value 

Primes axe used to indicate differentiation with respect to 
complete argument. Subscript notation is used to denote 
partial differentiation. 

METHOD OF APPROACH 

The configuration under consideration consists of a 
thin-walled unstiffened or ring-stiffened circular cylinder 
extending to infinity in the positive and negative £-direc- 
tions. (See fig. 1.) The cylinder is filled with a stationary 
fluid and is surrounded by air flowing in the positive 
r-direction at a Mach number M. The effects of midplane 
tensile stresses in both the circumferential and longitudinal 
directions and of a small amount of structural damping are 
taken into account. 




Figure 1 . — Infin itely long, thin-walled circular cylinder with air 
flow outside and stationary fluid inside. 


Por simplicity in the analysis, it is assumed that the de- 
formations of the cylinder walls can in most cases bo ade- 
quately described by Donnell’s differential equation. (See 
ref. 10.) It is kept in mind, however, that the validity of 
Donnell’s theory is limited to cases in which there are a 
large number of circumferential waves; where this condition 
is violated, a simplified version of FlOgge’s cylinder theory 
(see refs. 11 and 12) is employed to improve the accuracy 
of the results. 

The problem to be considered is the determination of those 
combinations of the parameters characterizing the cylinder 
and its environment that correspond to the boundary be- 
tween states of stable and unstable motion. For the pur- 
poses of this report, a system is considered stable if its 
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motion is either damped or purely sinusoidal; only timewise 
divergent motion is considered to be unstable. In order to 
determine the stability boundary itself, attention can be 
restricted to simple harmonic motion. However, such 
simple-harmonic-motion analyses often yield a multiplicity 
of boundaries, and it is necessary to derive equihbrium 
conditions for divergent oscillatory motions in order to 
determine the degree of instability in regions separated by 
the various boundaries and, thus, to identify the primary 
stability boundary. 

In line with the foregoing discussion, the method of ap- 
proach is first to derive the equilibrium conditions for si- 
nusoidal motion. These conditions are then extended to 
apply to divergent motion by means of analytic continua- 
tion. A Nyquist diagram technique is used to determine 
states of stability and stability boundaries. 

UNSTIFFENED CYLINDER 
DERIVATION OF EQUATIONS 

It is assumed herein that the cylinder wall may deform 
into any number of sinusoidal waves around its circumfer- 
ence and into sinusoidal waves of any wave length and 
constant amplitude along its length and, further, that the 
motion is simple harmonic in time; spacewise divergent 
motion (motion increasing in amplitude along the cylinder) 
is specifically excluded. Thus, the lateral deflection of the 
cylinder wall may be written 


Accordingly, the perturbation pressures exerted by a 
steady external or internal flow of Mach number Mi on a 
stationary cylinder deformed in the shape 


w=Ee(we _iKr i cos nd) 


( 2 ) 


have been determined in appendix A. The assumed defor- 
mation (eq. (1)) has the form of a wave traveling in the positive 
z-direction with velocity u>/v; hence, the flow of Mach number 
M outside the vibrating cylinder is equivalent to a flow of 
Mach number M — co/vc outside the stationary cylinder plus 
a flow of Mach number — w/rc* inside the stationary cylinder. 
Thus, by making the appropriate substitutions, the steady- 
flow results of appendix A are readily combined to give, for 
the outward air force exerted on a vibrating cylinder, 


Z(z,0,t)=flVB.e | ~lwe ^ ^ cos nd J (3) 

where 

7 = —pc?F(M—k, v, u)+p ( c, 2 i' , ( (j-Jc, v, nj (4) 
In equation (4), 


F(t,b,n)=? 


KJLt) 


5 zEY'{z) 



1 for &£> 0\ 

2 for &£<0/ 


(I5KD 

(HIM) 


( 5 ) 


w(x,d,t)~ Re(uie -, ’V"‘cos nO) 

=Re \_we v ” 1 cos nd] (n=2, 3, 4, ... ) (1) 

where w is the complex amplitude of the motion, v is the 
real wave number of the longitudinal waves, n is the num- 
ber of full waves around the circumference, and w is the 
frequency of oscillation. This assumed deflection shape 
will be the basis for the determination of the air forces 
exerted on the cylinder, the equihbrium condition, and, 
through it, the criteria for flutter. Note that n=0 and 
7? — 1 have been specifically excluded from consideration in 
this panel-flutter analysis because neither of these two 
motions involves panel action. The first value n=0 repre- 
sents pure dilation or contraction of the cross section; the 
second value n— 1 merely represents a rigid-body translation 
of the cross section. 

Air forces. — The air forces exerted on an infinite cylinder 
vibrating harmonically in still air have been reported in the 
literature. (See, for example, ref. 13.) Although some of 
these results can be applied to the present problem of de- 
termining the forces exerted on an oscillating cylinder in 
moving ah, it is convenient to derive this result directly. 
The unsteady-flow problem can be reduced to a steady-flow 
problem by means of a moving coordinate system in a 
manner similar to that of reference 7. 


z=#IV£ 2 -l 

In equations (5) and (6), J n is the Bessel function of the first 
kind, BY is the Hankel function of the first or second kind, 
and I n and K n are modified Bessel functions of the first and 
second kind, respectively. (See, for example, refs. 14 and 
15.) The primed quantity in the denominator of each func- 
tion is the derivative with respect to the entire argument. 
Note that a dimensionless frequency Jc=wjvc and a dimen- 
sionless wave number v=Rv have been introduced in equa- 
tion (4). 

Equilibrium condition. — Donnell’s equation for the equi- 
librium of thin, cylindrical shells (ref. 10) may be written in 
the modified form (see, for example, ref. 16) 


*«***{&> 
rf ^ "( g ) 

with 


(| 6 |< 1 ) 

(HIM) 


V (6) 


D 


V*w- 


1 12 ( 1 — M *) 
R* 






x ~~h 


£9 

R 2 


w e e'j+P,hw tt =l 

( 7 ) 
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where the subscripts x, 0, and t on w indicate differentiation 
with respect to these variables. The operators V 4 and V -4 
are defined by 



y 1 a» 


)' 


and 


V~ 4 V 4 =V 4 V~ 4 =1 


Substitution of w and l from equations (1) and (3) into 
equation (7) gives 


D 


(W- 


\ 


1 12(1— M 2 ) 


A 4 


p,b/j? — RiH = 0 


©' (^)'J 


+h 


as the condition for equilibrium of the motion of the cylinder 
wall in the presence of the moving air outside and the 
stationary fluid inside the cylinder. Written in terms of 
pertinent nondimensional parameters, this expression becomes 


X prr £&= 0 

M 


( 8 ) 


where 


0= 



12(1 At 2 ) r 


-Jm 


K 1 ^ 


- i n ~ 
'^iT- 2 ^ 


In the preceding definition, tr—a/E, and Cj—^E/p, is the 
speed of sound in the cylinder material. Note that, in 


equation (8), an additional term —i tk has been added in 

order to include qualitatively the effects of structural damp- 
ing of the Sezawa (viscous) type. The .damping coefficient 
£ is actually related to the parameters of the system in a 
rather complicated fashion but is always positive. A precise 
definition of e is not necessary because, in the following 
derivation, it is considered to be a vanishingly small, positive 
quantity. 


STABILITY DETERMINATION BY USE OF NYQUIST DIAGRAMS 

In order to determine the stability states of the structure 
whose behavior is governed by equation (8), it is expedient 
to examine the variation of the function L in the A-plano as 
k traverses, in the clockwise direction, a curve enclosing that 
portion of the & -plane corresponding to divergent motion of 
the cylinder (the lower half-plane for positive values of v or 
the upper half-plane for negative values of v). This, of 
course, is the well-known Nyquist diagram. (See, for 
example, ref. 7.) The number of resulting clockwise encircle- 
ments of the origin in the A-plane is equal to the number of 
zeros of L minus the number of poles of L enclosed by the 
curve traversed in the &-plane. Poles of L do not occur in 
the unstable half-plane since the functions comprising L are 
analytic everywhere in that region. (It should be noted 
that, in all instances, the real k ~ axis is excluded from the 
unstable half-plane.) Therefore, the Nyquist diagram gives 
directly the number of roots of equation (8) for which the 
frequency of oscillation has a negative imaginary part; in 
other words, it gives directly the number of modes of un- 
stable motion. 

With regard to the practical application of the Nyquist 
diagrams in the present case, it should be pointed out that 
actual computation of successive values of L is unnecessary; 
that is, through careful examination of the nature of the 
functions comprising L, it is possible to construct approxi- 
mate diagrams correct in all essential features and to see 
clearly the conditions under which critical changes (corre- 
sponding to flutter boundaries) occur. Further, it is only 
necessary to consider positive values of v because, if v is 
replaced by its negative and k by its conjugate, then L is 
replaced by its conjugate ; therefore, no new flutter boundaries 
would result from the negative values of v. 

ANALYSIS OF EMPTY CYLINDER 

Flutter criteria. — If the cylinder is assumed to be empty 
(j>t= 0), equation (8) reduces to 

L=— £ (Jt-Sl^-FiM-kyvri-i efc= 0 (9) 

p n |r| 


EXTENSION OF EQUILIBRIUM EQUATION TO DIVERGENT MOTION 

The flutter equation (eq. (8)) has been derived for sinu- 
soidal motion. As has been pointed out, an extension to 
divergent motion (of the form of eq. (1) with co a complex 
number having a negative imaginary part) is necessary. 
This requires the analytic continuation of equation (8) into 
the complex t-plane. The only terms in equation (8) for 
which this continuation is not trivial are the air-force func- 
tions F and F t . The analytic continuations of these func- 
tions are presented and discussed in appendix B. These 
functions are analytic throughout the half-plane correspond- 
ing to divergent motion and approach the values given by 
equations (5) and (6) as k becomes real. That the analytic 
extension of equation (8) into the lower half of the co-plane 
is indeed the equilibrium equation for divergent motion 
of the cylinder can be rigorously shown by the application 
of Fourier transform analysis. 


Let the damping coefficient e be very small but positive. 
With this restriction, the functions comprising L (for v]>0) 
vary with k in the manner shown in figure 2. Note that the 
inclusion of the infinitesimal damping is influential only in 
the range M—l<.k<^M+l where Im(— F)=0; elsewhere, 
Im(— F) is large in comparison with the damping term. In 
the range M— 1<C M+ 1 , the damping makes Im(Z) 

negative for positive k and positive for negative k. Hence, 
for subsonic flow (fig. 2(a)), the. damping insures that 
Im(Z)=0 only at k=0; for supersonic flow (fig. 2(b)), it per- 
mits Tm(Z)=0 only at k=M—l (in the limit as £->0). 
The significance of this will be made apparent by considering 
the variation of L as k traverses the path shown in figure 3 
(with the other parameters held constant). 

First, consider subsonic flow. (See fig. 2(a).) Two 
possible resulting variations of L are shown in the following 





(a) Subsonic flow (M < 1) . 
(b) Supersonic flow (M>1). 


Figure 2. — Typical variations with k of functions comprising L. 

(See eq. (9).) 


k- plane 


Figure 3. — Path enclosing lower half of fc-plane (corresponding to 
unstable motion for ?>0). 
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Nyquist diagram (where the full circle at infinity corresponds 
to the infinite semicircle in the Zr-plane) : 



For different values of the parameters, paths s imil ar to I 
or II may be traced. For path I, there is no encirclement 
of the origin and the cylinder is stable; for path Ft, on the 
other hand, one encirclement occurs and the cylinder js 
unstable. Since, by virtue of the damping term, Im(£) 
must pass through zero at &=0, it is apparent that the 
boundary between these two conditions is a static (diverg- 
ence) boundary defined throughout the subsonic range by 

Re (i)*_ 0 =0 = — ~ ~ Q 2 —F(M, v, n ) (10) 

Equation (10) may be put into the form 

< u > 

where 


_ 12(1 -M 2 ) 

“ 

1 

, - , n 3 - ~ 


(’■*) 

2 

5* 

K) 



and 

p_ 12(1— p 2 ) F(M,v,n) 

are positive numbers. Equation (11) has the one real root 

Mf+v « 

Thus, for selected values of the other parameters, the 
thickness ratio corresponding to the stability boundary in 




400194— C8 32 
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the subsonic range may be computed directly. This is the 
bounding value of h/R above which the cylinder is stable. 
For certain ranges of the parameters, namely when the 
Mach number is only slightly less than 1 or when the 
wave length of flutter is very large, the approximation 
M 2 

F(M, v, n ) « is permissible. When this approximation 

'TV 

is valid, equation (10) may be solved directly for M to yield 

(i3) 

Equation (13) gives the approximate hounding value of M 
below which the cylinder is stable. 

For supersonic flow (fig. 2(b)), the possible Nyquist 
diagrams resulting when k traverses the path of figure 3 
appear as shown in sketch 2. 



Since, now, the sign of the imaginary part of L always 
changes at k=M— 1, the criterion for a boundary is seen 
to be 

Re (L)*_ a r-!=0=^ ~ [(M— l) 2 — Q 2 ] — F(l, v, n) (14) 

Thus, the instability is dynamic, and, since ^ = ~’ the 
flutter mode is a traveling wave whose propagation velocity 
— is the velocity of external flow minus the velocity of sound. 

The solution of equation (14) for h/R can be carried out 
to yield a result in the same form as equation (12) with 
appropriate redefinition of A and B. In the supersonic 
case, however, the E-function is independent of M and is 
actually equal to —1/a; hence, equation (14) may be solved 


directly for the Mach number of flutter to give the convenient 
form 


M= lfl- 



dS) 


Determination of critical stability boundary. — In comput- 
ing stability boundaries, a factor which must bo kept in 
mind is the admissibility of all longitudinal wave lengths 
and of all integral values of n greater than or equal to 2. 
The critical boundary must necessarily correspond to the 
values of v and n for which the cylinder is most prone to 
flutter. This critical boundary can be obtained either by 
expressing the bounding value of h/R as a function of the 
rest of the parameters (as in eq. (12)) and maximizing this 
value of h/R with respect to both v and n or by finding the 
bounding value of M (as in eqs. (13) and (15)) and mini- 
mizing with respect to v and n. Of the two alternatives, the 
latter is more easily accomplished, analytical minimization 
of M being possible. For subsonic speeds, however, the 
expression for M (eq. (13)) is only approximate and, in some 
instances, it becomes necessary to maximize h/R as given 
by equation (12) . This is most readily achieved by graphical 
methods. 

Minimization of M can be performed by first minimizing 
n 2 with respect to v and then minimizing the resulting M 
with respect to n. The minimization with respect to v 
is particularly simple when there are no imposed midplalie 
stresses (<r x ~<rt—0); it can be performed with respect to 

the quantity 0+^0 • The result is 


tiff 


V12(1-m) 2 


at 


K) 


S Vl2(l— M 2 ) 

A 

R 


(16) 


(17) 


However, the quantity 0+S0 is itself a minimum at v=n; 
hence, equation (17) can be satisfied only when 


R 


(18) 


and, if this condition is violated, fl 2 can never achieve the 
value given by equation (16). The minimum value it can 
achieve is 


S2 


2 . 

min 


(c. Y 

IX®’- il 

W 

_12(1 — /j. 2 ) 1 4 n 2 _ 


(19) 


If the result given by either equation (16) or equation (19) 
is substituted into equation (15), the minimum value of M 
is seen to correspond to the smallest admissible value of n 
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(that is, 71=2) ; hence, for supersonic flow, the critical boundary is given by 




=1+ 


J Vl2(l- 

-p 1 ) 2 El — 
p R 



KDVl 

1 

\c) 

12(1 — p 2 ) 1 16 

2 p,h 

i 


p R 


(h ^ 712 ( 1 -jM 

VlT'- 16 ) 


> 


V120-M 2 ) 
16 


! ) 


( 20 ) 


Similarly, for subsonic flow, if the approximate equation, 
equation (13), is valid, the substitution of equation (16) 
into equation (13) yields the critical boundary 


M a 




(!)■ 


) V 12(1 — /t 2 ) 


, (1 -Me*) 


4 A 

* R 


visa-M 2 ) 


«a (2i) 


Numerical example. — For illustrative purposes, the critical 
stability boundary (n=2) has been computed for an empty, 
unstressed aluminum cylinder at sea level; additional curves 
corresponding to n=3, 4, 5, and 10 have been obtained for 
comparison. These results are shown by the solid curves 
in figure 4. Portions of the curves in the subsonic range 
were obtained by graphically maximizing h/R with respect 
to v as previously mentioned. It is interesting to examine 
the wave lengths of the flutter modes corresponding to these 
boundaries. In the range where equation (18) is satisfied, 
solution of equation (17) yields two different wave lengths 



Figure 4. — Stability boundaries for empty, infinitely long, unstiffened 
aluminum cylinder at sea level with no applied membrane stress. 


for the same critical thickness ratio. The larger one of 
these two wave lengths is shown in figure 5 for supersonic 
Mach numbers in the form of a plot of the aspect ratio of 
the flutter mode nfv (the ratio of longitudinal to circumferen- 
tial wave length) against Mach number for n= 2, 3, 4, and 5. 
The aspect ratios associated with the smaller wave lengths 
are not plotted; they are merely the reciprocals of the ones 
shown. At the higher Mach numbers, the critical value of 
h/R becomes large enough to cause equation (18) to be 
violated; when this happens, the two flutter modes coalesce 
to give a single flutter aspect ratio of unity as shown by the 
horizontal cutoff line in figure 5. 

As has been pointed out, Donnell’s theory is somewhat 
inaccurate for small values of n. (See, for example, ref. 12.) 
In order to obtain some idea of the magnitude of the resulting 
error, a new 12-function based on a simplified version of 
Fltigge’s cylindrical-shell theory (see, for example, refs. 11 
and 12) has been derived and minimized with respect to v; 



Figure 5. — Panel aspect ratio of critical flutter mode for empty, 
infini tely long, unstiffened aluminum cylinder at sea level with no 
applied membrane stress. 
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the details of this derivation are presented in appendix C. 
Computations in the supersonic range with the new (Flugge) 
O-function result in the dashed curves shewn in figures 4 
and 5. As in Donnell’s theory, the Flugge theory yields, 
for the lower thickness ratios, two values of wave length 
for which Q (and, hence, the Mach number) is a minimum. 
In contrast with Donnell’s theory, however, the two wave 
lengths are not equally critical; the larger one always yields 
a lower Mach number. For this reason, only the higher 
aspect ratio and its associated stability boundary are plotted 
in the figures. The critical boundary {n— 2) is found to 
require thickness ratios approximately 30 percent higher 
than those predicted by Donnell’s theory. The two curves 
for 71=3 still differ by 10 to 15 percent; but, for ti= 4, the 
two theories agree very well. For the sake of clarity, the 
stability boundary associated with the lower aspect ratio 
has not been shown in figure 4. It should be remarked, 
however, that this boundary agrees very well with the bound- 
ary given by Donnell’s theory, even for ti= 2. The prac- 
tical implications of the results shown in figures 4 and 5 are 
discussed subsequently. 


ANALYSIS OP FLUID-PILLED CYLINDER 


When the cylinder is assumed to contain a fluid, the addi- 
tional term F ( (^- k,v,n^ in equation (8) must be 


included. A typical plot of this function for real values of 
k is shown in figure 6. Note that the force is always real 
and becomes infinite at the resonant frequencies 



0 = 1 , 2 , 3 ,. . .) 


( 22 ) 


where the s/s are the zeros of J n ' (z). 

Since F { is always real for real values of k, the imaginary 
part of L can again be zero only at k=0 in the subsonic 
range and at k—M—1 in the supersonic range. Stability- 



boundaries can consequently occur only at these frequencies. 
In the subsonic range, the inclusion of F ( has no effect upon 
the stability boundaries since F t = 0 at k~0. In the super- 
sonic range, however, the inclusion of F t has a considerable 
effect on the stability boundaries, not only because of tho 
additional force at k—M—1 but also because of the reso- 
nances of F t . 

The condition for the stability boundary at supersonic 
speeds is 

Re(L) t . if _ 1 =0=^ | [(M-l)*-fl s ] + 

(23) 

This equation can be solved for h/R in the same form as 
equation (12) ; that is, 

h_ (B , [W~A?Y, (B IW 7T 3 Y S 

R V2~*V 4 + 27/ + \2 V 4 ~^27/ ^ 


where now 


12(1— n 2 ) 


0+0 LK’+S 


_ , 71?- (M-l)H 

2 -f7iT=j 


(?)‘J 


and 

B= 


I" l+Zi Y Ft (- (M— l),?,7iY] 

I I 7i^ Pj /c,^ \_n p \c / \Ci J _| 


The variation of h/R ivith Mach number for particular 
values of the other parameters is shown in sketch 3. Note 



Sketch 3. 

that h/R becomes infinite when the Mach number is equal 
tol+i,(i=l, 2, 3,. . . ). (See eq. (22).) This arises from 
the infiniteness of F t and, consequently, of B at those Mach 
numbers. Note also that the footpoints of the secondary 
stability boundaries occur when B=0. The numbers within 
the regions separated by the solid lines indicate the degree of 
instability as determined by use of Nyquist diagrams. 
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Since tlie value of h/R for the empty cylinder is given by 
equation (24) with F t equal to zero, and since F t is positive 
for values of M less than 1+Iq, it can be seen that the fluid 
inside the cylinder has a destabilizing effect. 

In order to find the critical stability boundaries, all values 
of v and n must be considered. Of particular interest in this 
connection is the fact that, for very large values of v, the 
infinities shown in the sketch all approach a value of M equal 

A. 

to l-( — . (See eq. (22).) Thus, the critical stability bound- 
c 

aries would appear as shown in sketch 4. For Mach numbers 



M 

Sketch 4. 


greater than 1 , the cylinder is unstable to an infinite de- 

c 


gree. If the fluid inside the cylinder is ah- at the same tem- 
perature as the surrounding air, this .limiting Mach number 
would be equal to 2. If the cylinder contains a relatively in- 
compressible fluid, however, this Mach number could be very 
high. In any event, the result is somewhat anomalous; it is 
probably caused by the use of linearized potential-flow 
theory and, undoubtedly, would not occur for real fluids. 


SOME REMARKS ABOUT THE SOLUTION FOR AN UNSTIFFENED CYLINDER 


On damping. — The stability criteria for the infini tely long, 
unstiffoned cylinder which were derived in the preceding sec- 
tions were obtained by including the effects of structural 
damping and then taking the limit as the damping ap- 
proached zero (e-H>0). This procedure was followed because 
different criteria are obtained when damping is not considered 
(esO). This important fact is illustrated in appendix D 
where the stability criteria for the empty cylinder with zero 
damping are discussed. Since structures always exhibit some 
damping, it is apparent that stability criteria obtained by 
taking the limit as e approaches zero axe more realistic than 
those obtained by setting e identically equal to zero. It is 
interesting to note that, in this case, the addition of damping 
makes the structure more prone to flutter. This somewhat 
surprising result may be explained by the fact that a damping 
force, even though in itself dissipative, can cause phase 
changes in such a manner as to allow the moving outside air 
to feed more energy into the structure; the result is a net 
energy gain. 


On the effects of midplane stresses. — The qualitative 
effects on flutter of the midplane tension stresses <r, and tr* 
(arising, say, from a static-pressure differential across the 
cylinder wall) may be determined by considering the defini- 
tion of- the wave speed Q yielded by the Donnell theory 
(immediately following eq. (8)). As expected, both cr x and 
tr e contribute to increased stability since they increase G and, 
hence, Mach number (see eq. (15)). Perhaps most signifi- 
cant, however, is the influence of midplane stress on the pre- 
dicted wave length of the flutter mode. Note that, while 
the presence of <s z stresses increases G uniformly for all wave 
lengths, the increase due to is proportional to the square 
of the aspect ratio n/v. Thus, the turn minima of G are no 
longer equally critical; the smaller minimum (and, hence, the 
most critical boundary) is seen to correspond to the shorter 
longitudinal waves. This is opposite to the result yielded by 
the more accurate Flilgge theory for the cylinder without 
midplane stress; thus, it would appear that addition of 
enough circumferential stress to the Flugge theory might 
produce a shift of the critical flutter mode from one with 
very long waves to one with very short waves. This possi- 
bility is important since results for an infinitely long cylinder 
would be applicable to a finite cylinder only if the wave 
lengths of the flutter modes were small in comparison with 
the length of the finite structure. 

Thus, in spite of the results yielded by the analyses of the 
unstressed cylinder pertaining to the wave length of flutter, 
it must not be concluded that analysis of an infinitely long 
cylinder has no application to cylinders of finite length. In 
fact, it may well be that, even for the unstressed cylinder, 
the minimum G corresponding to the predicted most critical 
flutter mode (with long waves) should be arbitrarily ignored 
in favor of the secondary minimum. Then, in view of the 
aforementioned agreement between Flugge’s and Donnell’s 
theories at this minimum, the Donnell theory would be ade- 
quate for all n(^2), the critical boundary being specified 
by the solid curve for n = 2 in figure 4. The aspect ratio 
of the flutter mode, in this case, would be the reciprocal of 
the value given in figure 5. 

On the effect of longitudinal stringers. — Although the 
analysis has been carried out for the case of an unstiffened 
cylinder, the flutter criteria may also be applied to the case 
of a cylinder with essentially rigid, longitudinal stiffeners. 
These longitudinal stiffeners would have the effect of raising 
the minimum value of n and, hence, of decreasing the criti- 
cal thickness ratio of flutter. For example, the curve for 
20 stringers (n=10) has been plotted in figure 4. 

This completes the discussion of unstiffened cylinders; an 
analysis of ring-stiffened cylinders is presented in the next 
section. 

RING-STIFFENED CYLINDER 

The ring-stiffened cylinder consists of the unstiffened cyl- 
inder with added, rigid ring stiffeners which prevent radial 
deflection at. the locations x=±ja (j= 0, 1, 2, . . .). (See 
fig. 1(b).) The stiffeners are assumed not to interfere with 
the flow of air outside or of fluid inside the cylinder. The 
analysis proceeds along the lines of that in reference 7. 
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DEBIT ATION OF EQUATIONS 

If the assumptions are made that the cylinder wall may 
deform into any number of sinusoidal waves around its 
circumference and into any shape periodic over 2 bay lengths 
in the longitudinal direction, the deformation may be 
written, with complete generality, as 


w(x, 6, t) 


=Ke( S 

\m— — c 


-i 


Omfi. 0 e tut cos nd 


) 


(25) 


provided that the coefficients a m satisfy the constraining 
relations 


S <u=o 

co 

S (— l)*a»=0 


(26) 


The functions F and F { are given by equations (5) and (6). 
They are repeated here for convenience. 




(lil<D 

_ t2 3?(z) 

/• f 1 for 6£>0\ 
\ \2for&£<[0/ 

<l€l>l) 



(ISKD 



(UI>D 


where 




These conditions correspond to zero deflection at the ring 
locations. They may also be written as 


S Ojn =0 


(m odd) 




S> =0 


(in acn) 


(27) 


Only one circumferential term is included in equation (25) 
because there is no structural or aerodynamic coupling be- 
tween the various cosine terms. Also, in this case, n—\ is 
admitted. The assumption of periodicity over 2 bay lengths 
is made because it permits a considerable degree of generality 
without overcomplicating the analysis. It is believed that 
the critical flutter mode would be of this type. 

Air forces. — In view of the linearity of the aerodynamic 
theory, the air forces exerted on a cylinder executing the 
motion given by equation (25) may be determined by 
separately considering each term of the s umma tion over m 
and superposing the results. Hence, the aerodynamic load- 
ing on a cylinder deforming in the shape of the traveling 
wave 

w m (x,e,t)—Re [*, 

is sought. 

If equation (28) is compared with equation (1), it becomes 
apparent that the air force is given by equations (3) to (6) 
when the substitutions w=a- m and v=mv/a are made. Thus, 
the total air force may be written immediately as 



2 / to -{'20. \ 

i(z,9,t)=l?^Re 0 e tat cos nBf (29) 

where, for m^O, 

l„= — n) 

(30) 


zHJIVFh 


In equation (30), &=coa/*-c and a= Rir/a. The quantity a 
is the ratio of one-half the circumference to the distance 
between the rings. For the special case of m= 0, a limiting 
process gives 

lo=—pc 1 G(k,a,n)+p i c i a Q i (j- k,ot,n) (33) 


In equation (33), 


G(Z,ct,ri)=? 


gyHti) 



1 for £<0\ 

2 for £>0/ 


(34) 


<?<(£,«,»)= £ 2 


JMi 1) 


(36) 


Equilibrium conditions. — By virtue of the assumed longi- 
tudinal spatial periodicity of the deformation, satisfaction 
of equilibrium over any interval of length 2a assures satis- 
faction of equilibrium over the whole length of the cylinder. 
Thus, it is sufficient to write Donnell’s equation for the 
segment 0^|a:<2a as 


D 


12 /Vx^ V-hiw 




1 — h We,') 


+ 


p,hw„=l(x,d,t)+P 0 (e,t)8{x)+Pi(6,t)8(x— a) (36) 


where Pq and P 1 are the reaction forces exerted on the cylinder 
wall by the two ring stiffeners included in the interval and 
where S(x) is the Dirac delta function. 

Let 


P 0 (8, t)=Re (Pot iu ‘ cos nd) 
Pi{6,t)='Re(P 1 e 1 ‘‘‘ cos nd ) _ 


(37) 


If P 0 and Pi from equations (37), w from equation (25), and 
l from equation (29) are substituted into the equilibrium 
equation, equation (36), and if this equation is then mill- 

imrx 

tiplied by « " and integrated over the interval, there 
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results 


2a < 


D 


/mV , n* V , 12(1 — n 3 ) ( 

V a 2 ^i? 2 / 


mirV 

. O’ / 


(c¥+m 




Hence, for p m ^0, 

where, for rn^O, 
with 


|a n =Po+(-l) s ’Pi 
(m=0, ±1, ±2, ±3, 


(hn z 


Po+(-l)"Pi 

2-ir > c 1 pR 


•) (38) 
(39) 


a 


P* 


«"-7 B S’ »«.»)+< |=I a 




C P m / ^/mo: 71 V 12(1— 1 . - . / B V_~ 

Vl2(l— /t^-y W V n -W + rt y+^+^ mg/ / ^ 


Note that, as in the case of the mastiff ened cylinder, a term 
has been added to each of the p m ’s to approximate the effect 
of damping (e„>0). In the limit as m-»0, 


/3o=-^|(fio 2 -F)-^0) 2 ft (|fc,a,»)+0(fc,«*,«)+i^fc 
where 



Finally, substituting from equation (39) into the constrairiing 
relations, equations (27), gives the following conditions for 
the existence of the motion defined by equation (25) : 

(Po-Pi) £ x-=0 (40) 

T71- — co Pm 

(to odd) 

(Po+Pi) £ 4-=0 (41) 

TO— ■ co p w 

(to <f«n) 

Inspection of these relations indicates that simultaneous 
satisfaction of the equilibrium equation and the restraint 
conditions can be achieved in a nontrivial fashion in several 
independent ways. One possibility is that 


Po=Pi 


(which implies that a„=0 when m is even) and 



(to odd) 


(43) 


As has been noted, the derivation herein and the results 
achieved exactly parallel those obtained in reference 7 for 
the infinite, flat plate. As pointed out therein, still another 
nontrivial solution is found when the restriction p m =^0 
(which characterized the foregoing results) is removed. 
Specifically, if two or more p m ’s of the same type (m odd or 
even) vanish simultaneously, a flutter mode may exist. 

Equations (42) and (43) are conditions for the existence of 
motion of constant amplitude; equation (42) corresponds to 
motion which is identical in each bay, whereas equation (43) 
corresponds to motion having the same amplitude from bay 
to bay but with alternating direction. 


STABILITY BOUND ABIES 


As in the case of the unstiffened cylinder, the Nyquist 
criterion can be used to investigate the states of stability of 
the stiffened cylinder. In this section are given the results of 
a limited investigation which was chiefly concerned with the 
examination of the following three-term approximation to 
equation (42): 


N= b w+ w - — 0 

Po PS P-s 


(44) 


(which implies that a m == 0 when m is odd) and 

£ f = 0 (42) 

TO-— o> p m 

(to aeri) 

and another is that 

Po=~Pi 


Subsonic flow. — In applying the Nyquist diagram tech- 
nique to the function N, account must be taken of the poles 
of N at the zeros of and p_ 2 . (There are no zeros of p 0 .) 
The results of such an application with e-»0 show clearly 
that, for subsonic Mach numbers, the only possible instability 
is divergence; further, the three-term approximation leads 
to as many as three static stability boundaries as shown in 
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sketch 5.~~PTn this sketch, a typical variation of ( h/R) max 



with M (with other pertinent parameters fixed) is shown. 
The subscript max indicates that the thickness ratio is 
maximized with respect to n. The numbers indicate the 
number of unstable roots of equation (44) corresponding to 
each region. Note that the upper curve corresponds to the 
condition that two or more jS m ’s of the same type vanish 
simultaneously, whereas the other curves result from the 
condition (A0t=o=0. 

Additional curves would result from the addition of more 
terms to N and still more curves would result from the 
consideration of the odd solution, equation (43); further, 
examination of a four-term approximation (m=l, —1, 3, 
and —3) to equation (43) indicates that, for reasonable 
values of a, the most critical boundary for either the even 
or the odd solution must stem from the condition 
03i)t«o=03_i)t=o=O- (For very small values of a, the dis- 
tance between rings becomes large and values of m other than 
unity may be critical.) The position of this critical boundary 
has been computed for an unstressed aluminum cylinder at 
sea level for various representative values of the aspect-ratio 
parameter a. The results are shown in figure 7. These 
curves correspond to various values of n ranging from 5 up- 
ward as illustrated on the plot; since n is large, Donnell’s 
theory is sufficiently accurate. 

Note that the thickness ratios required for static stability 
at subsonic Mach numbers are extremely small; therefore, 
divergence at subsonic speeds is probably not a critical 
design factor. 

Supersonic flow. — Application of the Nyquist criteria for 
supersonic flow is not so readily accomplished without the 
performance of further computations. It is clear, however, 
that divergence boundaries extend into the supersonic range 
and that flutter boundaries arise which probably become 
more critical than the divergence boundaries. 



Mach number, M 


Figtjbb 7. — Critical divergence boundaries for infinitely long, 
unstressed, ring-stiffened aluminum cylinder at sea level. 

The definition of these flutter boundaries requires that 
solutions be obtained to equations (42) and (43) or to suitable 
approximations, such as equation (44). This, in turn, re- 
quires extensive computations, especially since the resulting 
thickness ratios must be maximized with respect to n. This 
maximization can probably only bo achieved laboriously by 
making each computation for several values of n. 

CONCLUDING REMARKS 

A preliminary theoretical investigation of tho aeroelastic 
stability of infinitely long, thin-walled unstiffonod and ring- 
stiffened circular cylinders has been conducted by using 
Donnell’s cylinder theory and linearized unsteady potential- 
flow theory. A limited study of the resulting stability 
criteria has yielded the following information. 

For unstiffened cylinders with vanishingly small structural 
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damping, the only possible instability at subsonic Mach 
numbers is static divergence; in supersonic flow, however, 
flutter is found to occur for sufficiently thin cylinders in the 
form of a traveling wave whose propagation velocity is the 
velocity of the external flow minus the speed of sound. For 
an empty, unstressed aluminum cylinder at sea level, the 
critical boundary is found to correspond to a mode of de- 
formation having only two waves around the circumference 
but a very large longitudinal wave length. For this case, 
the use of the more complicated cylinder theory of Flugge 
is necessary for accurate determination of the critical ratios 
of cylinder-wall thickness to radius. These are found to be 
relatively large. 

Consideration of the effects of midplane tension stresses 
not only shows them to be stabilizing but indicates that large 
circumferential stress shifts the critical flutter condition from 
the mode with very long longitudinal waves to one with very 
short longitudinal waves for which Donnell’s theory is 
adequate. The magnitude of the stabilizing influence of 
the midplane stresses has not been investigated. 

The addition of an internal fluid has a destabilizing effect 


on the unstiffened cylinder. In fact, the anomalous result 
is found that for Mach numbers greater than 1 plus the ratio 
of the speeds of sound in the fluids inside and outside the 
cylinder, no adjustment of the physical properties of the 
cylinder will render it stable. 

The presence of even the smallest amount of structural 
damping is found to be an important factor in analyses of 
infinitely long, unstiffened cylinders. 

For the ring-stiffened cylinder it is found that flutter is 
not possible at subsonic Mach numbers and that only very 
small thickness ratios are required to prevent divergence. 
Although both panel flutter and divergence are possible at 
supersonic Mach numbers, no numerical results have been 
obtained; extensive computations would be required for a 
complete determination of the stability boundaries in this 
range. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., January 11, 1956. 
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APPENDIX A 

STEADY LINEARIZED FLOW PAST A STATIONARY DEFORMED CYLINDER 


In this appendix are derived the perturbation pressures 
exerted by a steady external or internal flow of Macb 
number M x on an infinitely long thin-walled cylinder with 
the deformation 

w='Re(we~ t ' Zl cos nd) (Al) 

where n is a positive integer and the wave number v may 
be either positive or negative. The Mach number M x is 
positive for flow in the positive Zi-direction and negative 
for flow in the negative r 1 -direction. (See fig. 8.) 

In the cylindrical coordinate system (x u r, 6), the linearized 
potential-flow equation is 

(1 — Mi)<t> ' +0rr+~ ^r+r 2 (A2) 

IjIj y y* 

where <t> is the velocity potential and the subscripts x u r, 
and 6 indicate differentiation with respect to these param- 
eters. With the assumption that 

^(r 1 ,r,0)=E,e[/(r)e -f ' It cos nO] (A3) 

equation (A2) becomes 

( f +i/ r +[,W/- I)-^j /=0 (A4) 

which, for supersonic flow (jM!|>l), has the solution 

<(r)=AB? (H +BH™ (|v| ^M^ALr) (AS) 



Figure 8. — Cylindrical coordinates for analysis of flow past 
stationary, deformed cylinder. 


and, for subsonic flow (|Mi|<l), has the solution 

f(r)=AI a (\v\ -Jl-M?r)+BK k (\v\ -fi=M?r) (A6) 

AIR FLOW OUTSIDE CYLINDER 

For air flowing over the outside of the cylinder, tho 
formulation of the problem is completed by tho specification 
of the boundary condition 

<t> r (xi Jt,d)=MiCW Xl (A7) 

at the cylinder wall and the proper conditions at infinity. 
By use of equations (Al) and (A3), equation (A7) becomes 

' f r (R)=—iMiCyw (A8) 

The resulting pressure perturbation A p can be calculated 
from the linearized Bernoulli’s equation 

Apfa ,r,e)=— P MiC<t> I1 (A9) 

and is 

Ap(x x , r,d)=pM 1 cRe(ifve~ irZl cos nO) (AlO) 

Supersonic flow. — For supersonic flow, the Sommerfeld 
condition (that is, that there be no incoming disturbances 
from infinity) requires that, for large values of r, the velocity 
potential <t> be essentially a function of Xi~ fir for M i>0 and 
a function of Xi+/9r for (where Then, 

from the relation for <f> (eq. (A3)), it is apparent that, for 

I Mi . r 

large values of r, j(r) must behave like e . Substituting 
the asymptotic approximations (ref. 14) 

1 (All) 

into equation (A5) gives, for large values of r, 

^( r )~V^/Sr [ e_ ^ n+ ^ 2 Ae 1| '^4-e^ n+5 ^J?«“ <l ^ r ] ( Al2 ) 

Hence, the Sommerfeld condition requires that B—0 for 
Miv^> 0 and A — 0 for M x v<C0. Therefore, equation (Afi) 
becomes 

HVHAlr) (i={2teM|v<o) <A13) 
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If, now, A is evaluated by use of the boundary condition, 
equation (A8), and the resulting expression for / is substi- 
tuted into equation (AlO), the perturbation pressure on the 
cylinder wall is obtained as 


Ap(x l ,R,6)=pc i 'Re^Rv i Mi 3 we -1 ’* 1 cos -nflj 

/ . Cl for M!» v>0\ 

V 1 2 for M!v<<y 

where 

z=E\v\^M l i -l 


(A 14) 


Subsonic flow. — For subsonic flow over the outside of the 
cyhnder, the potential must remain finite at infinity; hence, 
equation (A6) reduces to 


/(r)=5^(|r|Vr=M 2 r) (Al5) 

If B is evaluated by using the boundary condition, equa- 
tion (A8), and the resulting expression for /is then substi- 
tuted into equation (AlO), the resulting perturbation pres- 
sure on the cylinder wall is 


Ap(Xi,R,d)=pc i 'Re RAM 3 we~ i,z ' cos nO (Al6) 

L (D J 


where t=R\v\->]l— Mf. 

FLUID FLOW INSIDE CYLINDER 


at the cylinder wall is 

/ (R) = — iMiCfVW (Al 7) 

and Bernoulli's equation may be written as 

Api(x h r,d) = pMiCi ~Re(ifve- t ' x i cos nO) (A18) 

where the subscript i refers to the properties of the fluid 
inside the cylinder. 

For either supersonic or subsonic flow, the velocity poten- 
tial must remain finite throughout the cylinder. If this 
requirement is to be met at r=0, the solutions given by 
equations (A5) and (A6) must reduce, respectively, to 

' f(r)=AJ x (|v|VW=Ir) (AID) 

and 

m=AI n (\ v Wl=M?r) (A20) 

With the constant A determined through the use of the 
boundary condition, equation (Al7), and the resulting 
expressions for / substituted into equation (A18), the per- 
turbation pressure for |Afi|)>l becomes 

Api(xi, R, 0)=p 1 c, i Re we~ irZl cos nflj (A21) 

and, for IMjKI, 

Ap t (x 1 ,R,d)=p i c/Jic j^7?F*Mi 2 we~ irZl cos ntfj (A22) 


It is assumed that the fluid inside the cylinder may be 
other than air and may exist under different conditions from 
the air outside the cylinder. Thus, the boundary condition 


These results hold for flow in either direction. Note that, for 
these air forces, the pressure is always in phase with the 
motion. 
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APPENDIX B 


AIR FORCES ON A CYLINDER IN DIVERGENT OSCILLATION 


The outward lift force on a cylinder wall executing the 
simple harmonic motion (« real) 

w(x,6,t)='Re(we~ irX e iut cos n8) (Bl) 

in the presence of moving air outside and a stationary fluid 
inside the cylinder has been determined. (See eqs. (3) to 
(6) .) The amplitude of this force is expressed by 


7= - p( ?F (M~, By, 7 ^+ Pi c t 2 F t (A Rv,n) ©2) 

In equation (B2) , 




(l*ld) 

- S^iz) - , 
zH<?'(z) \ 

f. r 1 for i?j<£>0\ 

* = \2 for R^COj 

(UI>D 

and 





(|{|<D 

JO/ 

II 


(lfl>D 


K(B3) 


r (B4) 


where 


r=i?MV i-i 5 


z = R \ v \^¥ — 1 


It is desired to extend these air-force functions to apply to 
divergent motion of the cylinder (u complex -with a negative 
imaginary part). This extension can be made by using 
analytic continuation; in fact, it can be shown by Fourier 
transform analysis that the necessary and sufficient condi- 
tions for the extensibility of the air-force functions to 
divergent motion are that the air-force functions possess 
analytic continuations which have no singularities anywhere 
in the unstable half-plane (Im(co)< 0) and reduce smoothly 
to equations (B3) and (B4) on the real axis. Singularities 
may possibly exist along the real axis, but the real axis is 
specifically excluded from the unstable half-plane. 

Consider the function F(£,Rv,n). Examination of the 
manner in which co appears in F (see eq. (B2)) shows that the 
unstable half of the f-plane is the upper half-plane if v is 
positive and the lower half-plane if v is negative. It can be 
verified that the desired analytic continuation of F into the 
proper half-plane is given as follows: 


F(£, Rv,n)= 


ggnCgj HVEg) 


(B5) 


where the desired branch of the multiple-valued function F 
is the one for which the cuts are as shown in sketch 6: 




v> 0 v<0 

Sketch 6. 


and 


F(V 2, Rv, n) — 


2H?{R\v\) 

R\ v \H?'(R\v\) 


R\v\nrm) 


(?> 0)' 

«0) 


(B6) 


As Im(£) approaches zero from the proper direction, tlio 
definitions given by equations (B3) result. Furthermore, 
no branch points or cuts of the function F occur in the 
unstable half of the f-plano, and, in addition, Nyquist 
diagram techniques can be employed to show that there are 
no poles of Fin this region. 

In a similar manner, the analytic continuation of F t can 
be expressed in the form 


F { (£,Rv,n)= 


?j,CsMvr=?) 

R\v\JT=? I n '{R\v\J[=?) 


(B7) 


No difficulties with branch points occur for F t ; the functions 
I n and I B ' are entire functions and the combination IJzIn 
is an even function. Therefore, F t is inherently a single- 
valued function of £. 

The air forces used in the analysis of the ring-stiffened 
cylinder can be similarly written for divergent motion. For 
m 7 * 0, it is only necessary to replace v by rmr/a in equations 
(B5) to (B7). For m= 0, it is necessary to specify the air 
forces separately. Thus, for divergent motion (seo eqs. 
(33) to (35)), 


where 


l 0 =—pc 2 O(k, a, 1c, a, nj 


n)= 




(B8) 

(B0) 
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(B10 > 

The desired branch of the multiple-valued function 0 is the 
one for which the cut extends from the origin to infinity 


along the positive imaginary axis and 


G(X,n,a)=z 




(Bll) 


APPENDIX C 

EQUILIBRIUM OF UNSTIFFENED CYLINDER— FLtlGGE’S EQUATION 

The simplifying assumptions used in deriving Donnell’s differential equation (eq. (7)) result in a loss of accuracy 
when the condition n 2 »l is violated; further, terms which become important when the longitudinal wave lengths 
become large (that is, when the behavior of a cylinder approaches that of a ring) have been omitted. 

The more complicated theory proposed by Fltigge (see, for example, ref. 11) is not characterized by either of these 
limitations. FlOgge’s equations may be written in the form of a single equation in w, which, for (h/R) 2 <^<^l 
reduces (see refs. 11 and 12) to 

D^V 4 w+-~ V~ 4 w xzSS +~ w em + 2iiRhu xxm +Qw zxxz «B+ 2( ~j^ ^ w xxlte$ +~ -+p,hw tt =l (Cl) 

where imposed midplane stresses have not been included. As in Donnell’s equation, equation (Cl) takes no account 
of inertia forces in the longitudinal and circumferential directions. 

Substitution of w and Z for the unstiffened cylinder (eqs. (1), (3), and (4)) into equation (Cl) gives 



Equation (C2) is precisely equation (8) ; hence, the flutter boundary in the supersonic range for the empty, unstressed cylinder 
is still defined by 


(C4) 



but with 9 given by the more complicated function, equation (C3). As before, the critical boundary corresponds to the 
minimum value of M as v and n are varied. Setting b9?jbv= 0 yields 


V°+Zti?P+\ 2» 4 +2(3— 2Dn?- 
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-Zn\n 2 -Vf?-n\n 2 -l) 2 =0 (C5) 


Careful examination of equation (C5) indicates that, for small values of h/R, there are two minimums of 9 2 with an intervening 
maximum. Computations have been made for an aluminum cylinder at sea level by letting n= 2, 3, and 4; both the minimums 
of 9 2 were chocked and, through equation (C4), the boundaries corresponding to the lower of these were determined. These 
results are shown by the dashed curves in figure 4; in addition, the corresponding wave lengths are given in figure 5. 
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APPENDIX D 


IMPORTANCE OF DAMPING IN UNSTIFFENED CYLINDER 


It is of interest to examine some of the results of an 
analysis of the empty, unstiffened cylinder which takes no 
account of damping (as opposed to the analysis in the body 
of this report which treats a cylinder in the limit as damping 
approaches zero). The equilibrium condition, equation 
(9), reduces to 


Le-A (fc*_ v,n)=0 (Dl) 

p K 

and its constituent functions vary as in figure 2 except that 
now Im(L) =0 throughout the range M — 

A typical Nyquist diagram (resulting when k traverses 
the path of fig. 3) is shown in sketch 7 which corresponds 



to J>he particular conditions ?>0, 0<M— 1<(&<\M, 
Re(A) t .af_i|>0, and Re(dA/di)j taJ f_ 1 <(0. (The infinitesimal, 
counterclockwise semicircles correspond to infinitesimal, 
counterclockwise semicircles traversed by k to exclude zeros 
of L on the real axis from the lower half-plane.) These 
particular conditions are illustrated because they show 
clearly the possible existence of a stable condition with 


R.e(X-)jt»Af_i>0, a circumstance found impossible with damp- 
ing present. The reason for the difference, of course, is 
that, for &>0, the damping term is a negative imaginary 
quantity and even the smallest amount of damping shifts 
all the portion of the diagram corresponding to the range 
M— l<A<(M-fT below the real axis. 

Note that, for the case of zero damping, satisfaction of 
equilibrium (X=0) with k real is not always sufficient for 
definition of the boundary; for the conditions in the cited 
example, the roots of L= 0 do not pass directly from tho 
upper into the lower half of the i-plane as the other param- 
eters are varied, but linger awhile on the real axis. In- 
stability occurs when a root leaves the real axis and enters 
the lower half of the -plane. For the cited example, this 
definition corresponds to the simultaneous solution of the 
equations 

1 = 0 } 



(D2) 


in the range M— 1<£<M. 

Another significant difference in the two analyses is that, 
when damping is ignored, the instability in the subsonic 
region is low-frequency flutter and no purely static instability 
of the cylinder is possible. Similar differences would bo 
found for the fluid-filled cylinder if the damping were 
ignored. The fact that the two assumptions, e->0 and 
6=0, produce such discontinuously different results testifies 
to the importance of tho presence of damping in tho un- 
stiffened cylinder. 
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